Charge Transport Scalings in Turbulent Electroconvection 
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We describe a local-power law scaling theory for the mean dimensionless electric current A'^m in 
turbulent electroconvection. The experimental system consists of a weakly conducting, submicron 
thick liquid crystal film supported in the annulus between concentric circular electrodes. It is driven 
into electroconvection by an applied voltage between its inner and outer edges. At sufficiently large 
voltage differences, the flow is unsteady and electric charge is turbulently transported between the 
electrodes. Our theoretical development, which closely parallels the Grossmann-Lohse model for 
turbulent thermal convection, predicts the local-power law Nu ~ F{r)TZ"'P^ . TZ and P are di- 
mensionless numbers that are similar to the Rayleigh and Prandtl numbers of thermal convection, 
respectively. The dimensionless function F{T), which is specifled by the model, describes the de- 
pendence of Nu on the aspect ratio F. We find that measurements of Nu are consistent with the 
theoretical model. 
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INTRODUCTION 

Turbulent Rayleigh-Benard convection (RBC), the 
paradigm for studies in convective turbulence, remains 
a fascinating unresolved puzzle 1]. For a century, theo- 
ries of turbulent RBC have focused on understanding the 
globally averaged heat transport through a layer of fluid. 
The development of models has been largely driven by 
improved experimental measurements of the heat cur- 
rent, which have time and again revealed unexplained 
discrepancies between experiment and theory. In re- 
sponse, mathematical models have become increasingly 
sophisticated during the last five decades 0. In recent 
years, ambitious experimental, theoretical and compu- 
tational projects have been undertaken and the study 
of turbulent RBC has been considerably reinvigorated 
^, 4, 5, 6,, 7 , 9, iO, H, il, 13, li, OSnlfi] . All this ac- 
tivity suggests some clear directions for future work. Two 
recent experiments highlighted the crucial role of the sys- 
tem shape and lateral extent 9, 16]. In this paper, we 
exploit the unique features of a system closely analogous 
to RBC, electrically driven convection in a thin annular 
film, to shed light on these features of turbulent convec- 
tion. 

Turbulent RBC is described phenomenologically in 
terms of several organizing structures which are exper- 
imentally observed. The convecting fluid has sharp ther- 
mal boundary layers at its top and bottom surfaces. 
Plumes grow erratically from these surfaces and spon- 
taneously organize into a noisy but coherent wind. This 
large scale circulation (LSC), or turbulent wind, advects 
the plumes so that hot thermals rise along a laterally 
bounding wall while cold plumes sink along the diamet- 
rically opposite wall. The interior of the cell away from 
the LSC is typically assumed to be well-mixed, isotropic 
and homogeneous. This picture is partially corroborated 
by measurements using containers with approximately 
equal height and breadth. The extent to which this phe- 
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FIG. 1: Schematics of the experiment: (a) top and (b) side. 

nomenological picture, particularly of the LSC, is geom- 
etry and shape dependent is currently the most pressing 
question in the study of turbulent RBC. 

The aspect ratio F, which is the ratio of the horizon- 
tal span to the vertical separation of an RBC apparatus, 
quantifies the geometry. It is important to perform ex- 
perimental studies with a wide range of F. These have 
been difficult to accomplish however, due to practical 
considerations. The requirement of strong forcing, which 
favors large vertical separation, is opposed by the con- 
straint that the apparatus remain a manageable labora- 
tory size, limiting its horizontal span. Consequently the 
majority of experiments have F w 1. Turbulent elec- 
troconvection, which is largely unencumbered by scale 
considerations, has emerged as a complementary exper- 
imental system for the study of convective turbulence 

Electroconvection driven by DC potentials in thin 
smectic liquid crystal films has been extensivel y st ud- 
ied, mostly in the weakly driven, laminar regime |l7lll8l 
El El EE m f23, 24, 25, 26, 2LM- The system con- 
sists of a smectic A liquid crystal film suspended in the 
annulus between the edges of concentric circular metal- 
lic electrodes as shown schematically in Fig. ^ An ap- 
plied electric potential difference between the inner and 
outer electrodes drives an electric current through the 
film. Surface charges accumulate on the two free sur- 
faces that separate the electrically conducting film from 
charge-free space. This inverted surface charge density 



is unstable to electric forcing in much the same manner 
as the inverted mass density distribution of RBC is un- 
stable to buoyancy forces. When the applied voltage V 
exceeds the critical voltage Vc, the fluid is organized into 
convection vortices as shown in Fig. The flow ad- 
vects electric charge between the electrodes constituting 
a convection electric current. At higher driving, the flow 
becomes unsteady while retaining the large scale struc- 
ture of convecting vortices, as shown in Fig. Eh. The 
turbulent electric charge transport is analogous to the 
heat flux in turbulent RBC. 

Other AC driven forms of electroconvection have been 
studied in bulk, three-dimensional liquid crystals [IM 113 
In this paper, we concern ourselves only with 
two-dimensional, surface charge driven convection under 
DC potentials. However, it seems clear that the turbulent 
regime in these other systems might be amenable to a 
similar sort of scaling analysis . 

The globally averaged heat transport in RBC is quan- 
tified in dimensionless form by dividing out the contri- 
bution due to molecular conduction. It is then referred 
to as the Nusselt number Nu. Almost every theoretical 
model over the last 5 decades has attempted as its cen- 
tral goal to describe the functional dependencies of Nu. 
It is generally accepted that Nu — Nu{Ra, Pr,T). Here 
Ra is the Rayleigh number which quantifies the thermal 
forcing, Pr is the Prandtl number which is the ratio of 
the fluid's kinematic viscosity to thermal diffusivity and 
r is the aspect ratio of the system. 

Early models for Nu used the laterally extended limit, 
r — > oo, while some recent theories assume the con- 
fined case of r « 1 0. A few years ago Grossmann 
and Lohse (GL) proposed a unifying theory for the Ra 
and Pr scaling of global variables like TV it [6j, albeit for 
the case of F « 1. The GL scheme begins by decom- 
posing thermal turbulence into two constituents: a ki- 
netic component and a thermal component. Then, us- 
ing a second phenomenological distinction between the 
boundary layer and the bulk, the GL theory estimates 
the kinetic and thermal dissipations in each region us- 
ing dimensional arguments. Finally, GL theory derives 
the dependence of Nu on Ra and Pr by balancing the 
exact formulations for the total dissipations versus the 
dominant contributions from the bulk and/or bound- 
ary layer. The various combinations of bulk/boundary 
layer, kinetic/thermal, and other considerations lead to 
the fragmentation of the Ra — Pr parameter space into 
ten regimes 0- Within each regime, GL theory predicts 
the local-power law scalings Nu ~ Ra'^Pr^ with regime- 
dependent exponents 7 and 6. Near regime boundaries, 
we expect crossover effects. Unlike previous scaling the- 
ories, GL theory predicts no purely power law scaling for 
Nu. Rather, Nu — f{Ra,Pr) can be approximated by 
regime-dependent combinations of local power laws for 
F « 1 systems. Recent precision experiments are bet- 
ter modeled by the GL cross-over scaling theory than by 




FIG. 2: (Color online) Qualitative visualization of electro- 
convective flow, characterized by a fllm of uneven thickness: 
At (a) V = 100 volts and (b) V = 250 volts. Data, however, 
were acquired only from uniformly thick films in which no 
flow is visible. 

simple power laws 0, E| ■ 

Deep similarities between the phenomenologies and 
mathematical descriptions of turbulent electroconvection 
and RBC make possible the development of a GL the- 
ory for electroconvection. In this paper, we exploit these 
parallels to describe the nondimcnsional electric current 
(also denoted Nu, the electric Nusselt number) in turbu- 
lent electroconvection. We find that Nu = F{T)f{TZ,V) 
where TZ and V are the electrical analogs of the Rayleigh 
and Prandtl numbers of conventional RBC. The func- 
tion f(Ti,V) can be approximated by local-power laws 
with the same regime dependent exponents as in the GL 
theory. We find that the theoretical predictions are con- 
sistent with electric current measurements. Unlike the 
GL model, where it is assumed that F « 1, the naturally 
periodic annular geometry of electroconvection permits 
the explicit calculation of the aspect ratio dependence 
F(r). The aspect ratio dependence of Nu is not a power 
law as is the case in the theory of Shraiman and Siggia 
[ssj i. but instead approaches a constant for large F. We 
find surprisingly broad agreement between the function 
F and measurements from both turbulent electroconvec- 
tion and RBC 
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This paper is organized as follows. We start by describ- 
ing the physical model for annular electroconvection in 
thin freely suspended liquid crystal films. We then de- 
rive rigorous results for the globally averaged kinetic and 
electric dissipations. Using two examples that are di- 
rectly relevant to the current experiments, we show the 
derivation of local power law scaling for Nu and for the 
Reynolds number Re of the large scale circulation. Im- 
portant features of the experimental system are then re- 
counted and we describe direct comparison between the 
scaling theory and experimental data. We finish with 
a brief conclusion where we discuss the complementary 
relationship between our work on turbulent electrocon- 
vection and conventional RBC. 



SCALING THEORY 
Background 

The basic equations for annular electroconvection are 
introduced in this section; for full details see Ref. 
A close parallel with RBC is evident, which allows us 
the develop a scaling theory, analogous to GL theory, in 
subsequent sections. 

The annular film is a two dimensional (2D) sheet 
that in a cylindrical coordinate system spans the region 
ri < r < To at z = 0. Here (r,,) are the radii of the 
inner (outer) electrodes. On the plane z = 0, the in- 
ner (outer) electrodes occupy the region < r < Vi (r,, < 
r < oo). The film is a Newtonian, incompressible, elec- 
trically conducting fluid with 2D density p, molecular 
viscosity rj, and electric conductivity a. We write the 
equations governing the fluid and charge in the annular 
film using 2D differential operators, field variables and 
material parameters. The fiuid fiow is described by the 
Navier-Stokes equation with an electric body force: 



u 



Vu) = -VP 



(1) 



Here the fluid velocity u = {u,v), P is the pressure, q the 
surface charge density, E is the electric field and gE is 
the electric body force. The conservation of charge leads 
to a continuity equation: 



dtq + \7-J = 0, J = CTE + gu. 



(2) 



The current density J is composed of the usual Ohmic or 
conduction density crE and the convective current density 
qu. The electric field and charge density are not indepen- 
dent variables and have to satisfy Maxwell's equations. 
Magnetic effects are negligible. In the region \z\ ^ i.e. 
above and below the film and electrodes, space is free of 



charges and so Laplace's equation for the electric poten- 
tial holds: 



V3V3 



0, V3 



v2 + 5. 



(3) 



In the above, the subscript 3 identifies the potential and 
gradient operator are defined in three dimensions. The 
potential "03 can determined in the upper half space by 
solving Eqn. 01 subject to boundary conditions at z = 
and at infinity. The surface charge density q on the film 
due to the discontinuity in the electric field normal to the 
film is given by 



-2eo9^V3U=o+ 



eo9z'03|z=o- 
r,; < r < To 



(4) 



In the above eo is the permittivity of free space. Equa- 
tions QEl model the electroconvection system. The equa- 
tions are subject to rigid boundary conditions on the fiuid 
velocity: u = at r = and r — Tq and to Dirichlet 
boundary conditions for the electric potential at z = 
and at infinity. The applied potential is V volts for r < Vi 
and volts for r > Tq and at infinity. On the film, the 
electric potential is determined by satisfying the current 
density J and the boundary conditions at r = r^, Tq. 
Using the Dirichlet Green function G{r, 6, z; r' , 0' , z'), we 
can formally solve Eqns. |31and0]for the surface charge 
density: 



ZTT 



da' il;3{r',e',z' = 0) 
g[ilj], ip = il;3{r,0,z 



dzdz' 
0). 



z=0+ 



(5) 



The integral is over the bounding surface, in this case 
the z — plane and surface at infinity. Eqn.[3essentially 
defines q in terms of a functional g whose argument is 
the electric potential in the plane z — 0. On this plane 
the electric field E and potential tp are related though 
E = —Vip. Using this relation, the definition of g, de- 
noting the kinematic pressure field p = P/p and kine- 
matic viscosity v = ri/p, the four Eons. ITEI reduce to the 
following pair: 



5tU-^u-Vu = -Vp-hJ/V^u- — V?A, (6) 

P 

StV + u-W = (7) 

o^^g 

Written in this way, the Eqns. and [7| bear striking 
similarity to the Boussinesq equations for turbulent RBC 
with the scalar temperature and electric potential fields 
assuming similar roles. The solutions of Eqns. and [7| 
are subject to the usual no-slip boundary conditions on 
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u and the applied electric potential boundary conditions 
on ^p. Three dimensionless parameters describe the state 
of the system. TZ, the analog of the Rayleigh number in 
RBC, is the control or external driving parameter, and 
is proportional to the square of the applied voltage. It is 
given by 



n 



(8) 



V, the analog of the Prandtl number in RBC, is the ratio 
of the charge to viscous relaxation time scales in the film: 



V 



£0V 



The geometry is uniquely characterized by the radius ra- 
tio a = Til To. However, in order to make comparison 
to RBC, it is more appropriate to describe the geom- 
etry in terms of the aspect ratio F, which is the ratio 
of the horizontal or lateral dimension to the vertical or 
transverse dimension. Since the lateral dimension is am- 
biguous in the annular geometry, F is not uniquely deter- 
mined. However, F can be consistently taken to be the 
ratio of the circumference of the film measured at the 
inner electrode to the film width, so that 



F = 



1 -a 



(10) 



Two possible alternative definitions of the aspect ratio 
use the mid-radius circumference or the outer electrode 
circumference as the lateral dimension, leading to F„j and 
Fo, respectively. The two aspect ratios so defined are re- 
lated to F defined in Eqn. ^jas F = F„i — -k — V o — Stt. 
Since the various definitions of the aspect ratio are very 
similar, we have chosen to use the form in Eqn.^l be- 
cause in that case < F < oo, as is true for conventional 
RBC systems. The alternative definitions F^, and Fq 
have TT and I-k as their lower bounds. In any case, direct 
comparisons between rectangular and annular systems 
converge in the a ^ 1 or F — > cxd limits. 



Global averages 

In this section we derive exact expressions for various 
globally averaged quantities, which, according to the GL 
procedure, will be balanced against scaling estimates, as 
discussed in the following section. 

The convective contribution to the electric current is 
determined by dividing out the conduction current from 
the total current. This dimensionless ratio is the Nusselt 
number Nu. The net electric current is radial between 
the inner and outer electrodes with a net zero contribu- 
tion from azimuthal currents. Thus at any radial position 
Ti < r < To, the integral over the azimuth of the radial 



component of the current density is equal to the total 
electric current. Hence, we find 



Nu = 



■J-rdl 



§ uq — adj-ip dl 



V dl 



-crdrtpc 



dl' 



(11) 



In the above, we have used the definition of J given in 
Eqn. 121 The contour dl is a circle at radius r, and f is 
a radially outward unit vector. Note that u is the radial 
component of the fluid velocity u and ipcond is the electric 
potential for the conductive state. On the inner electrode 
(r < Ti) ipcond = V while on the outer electrode (r > Tq) 
V'cond = 0. On the film where n < r < Vo, ipcond is given 



(9) by 



■0c 



V , r 
In — 

in a To 



(12) 



Evaluating the terms in Eqn. ^2 using Eqn. ^1 we find 
the relation 



Nu : 



J^^{uq — adrtp) rd9 
27r(TF/ln(l/a) 



(13) 



We use the angular brackets (• • •) to denote averages 
over the fluid volume (actually its area), we write 



(...) 



/ • • • rdrde _ Jo"" • • • rdrde 



J rdrdO 



TT{rl - r2) 



(14) 



As in the GL theory for RBC, our charge transport scal- 
ing theory begins with the kinetic and electric dissipation 
rates 



(15) 



We denote the averages of the above dissipations over the 
fluid volume as 



eu = (eu(r,0,t)) = (^^(Vu)^) 
= (e^(r,0,t)) = (a(V7A)2) 



(16) 



The rather complex and nonlocal functional relationship 
g[ip] between q and ip given in Eqn. |S1 makes it very dif- 
ficult to calculate the globally averaged kinetic and elec- 
tric dissipations. However, by splitting the integrand in 
Eqn. 13 into local and nonlocal parts, g can be expanded 
as 



5M = 



-ijj + nonlocal terms . 



(17) 



To a first approximation the functional g[ip] can be ap- 
proximated as a linear function g{ip) so that the surface 
charge density and the electric potential on the film are 
related locally. Then d^g is constant and Eqn.[71is iden- 
tical to the heat equation in RBC, with electric potential 
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in place of temperature. The critical parameters and the 
critical mode numbers at the onset of electroconvection 
have been successfully captured by this local approxima- 
tion, which is described in detail in Ref. . We assume 
that g (X ip continues to hold at much higher forcing 
and develop the GL theory for turbulent electroconvec- 
tion within the local approximation, which makes Eqns.|5| 
and 13 very similar to the Boussinesq equations for ther- 
mal convection. The electric potential in electroconvec- 
tion almost precisely corresponds to the temperature in 
RBC. 

The velocity and electric potential fields share the az- 
imuthal periodicity of the annulus and thus permit the 
calculation of the globally-averaged kinetic and electric 
dissipations of Eqn. 1161 Assuming time-stationarity for 
the spatial averaging, we find from Eqns.ElandlZlthc fol- 
lowing relations for the kinetic dissipation eu = (i^(Vu)^) 
and electric dissipation = (fT(VV')^): 



t/37^p-^(Nu- 1) 
ln(l/a)(r2 - rf){ro - 
2ct1/2Nu 



(18) 
(19) 



The above relations have interesting similarities and 
differences to those for the corresponding quantities in 
RBC, which are given in Refs. j^l^jH^- The differences 
are entirely due to the annular geometry of the electro- 
convection system. In the narrow gap limit, in which the 
radius ratio a — > 1 while film width d = Tq — remains 
constant, Eons. O [THl and [T^ recover the familiar forms 
for RBC between parallel plates. 



Grossmann-Lohse scalings 

In this section, we make several assumptions about the 
spatial organization of the turbulent flow, in the same 
manner as in the scaling theory of Grossmann and Lohse 
(GL), as they explain in Section 2 of Ref. 0. In particu- 
lar, GL assumed that a turbulent wind, or LSC, compris- 
ing a single cell occupies the entire F ~ 1 RBC container. 
The wind is driven by plumes from the boundary layers 
and it in turn drives the interior or bulk. With the bound- 
ary layers and bulk conceptually distinguished, GL esti- 
mate the relative boundary and bulk dissipations. Here, 
we make analogous assumptions about turbulent electro- 
convection. As shown in Fig.|23, the turbulent flow con- 
sists of counter-rotating convection vortices around the 
annulus. The vortices, which are unsteady, have fluctu- 
ating boundaries that are defined by the averaged turbu- 
lent LSC. Each vortex is assumed to be roughly square 
with dimension — ri. Near the electrodes, well devel- 
oped viscous and electric boundary layers with respec- 
tive thicknesses Au and are assumed. Away from 



the electrodes, the vortex interior, or bulk, is taken to 
be well mixed. We do not account for the slight differ- 
ences between the boundary layer dimensions at the inner 
and outer electrodes due to the annular geometry. This 
asymmetry diminishes with increasing aspect ratio and 
is assumed to be always small. 

The total dissipations calculated in Eqns. ^1 and [T^ 
are decomposed into contributions from the boundary 
and bulk regions of the convection cells as follows 



hi . bulk 

^11 t,, I c„ 



hi I hulk 



(20) 



The contribution ej,' of the boundary layer kinetic dissi- 
pation is defined as 



(eu(r-,6l,i))fci = {v[VvLf)u, 



2tt rr-i + Au 



/o /, 



• rdrdO 



p27r pro 
JO Jr. 



■ ■ ■ rdrdO 



Arl - (ro - A„)2) 



(21) 



The electric dissipation in the bulk e^"'*^ is defined as 

- {f-^{r,9,t))buik ^ {(j[^i>f)buik , 



hulk 



(...) 



2-K rf^o — ^Ti 



rdrdO 



bulk 



7r((ro - A^)2 - (r, + A^)2) 



(22) 



The other two boundary and bulk dissipations are simi- 
larly defined. We assume that 



Nu = F(r)f(7e,7') 



(23) 



It is our purpose to determine the as yet unspecified func- 
tions F and f. We begin by estimating the bulk and 
boundary layer dissipations. In the following we ignore 
numerical factors of 0(1). 



Kinetic dissipations 

The turbulent wind or LSC sets the velocity scale U 
for both the boundary and bulk regions. A viscous or 
kinetic boundary layer, assumed to be laminar, scales as 
Au ~ (fo — ri)Re~-^/^. Here Re is the Reynolds number 
based on velocity U and film width Tq — r.^. Finally we 
assume that Au «C J'i, Au <^ and Au <IC — which 
greatly simplifies the averaging in Eons. ET1 and 1221 For 
the boundary layer, we estimate the kinetic dissipation 
by 



A,: 



Au / To ^ Ti 



{ro - n)- 



(24) 
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The LSC stirs the interior and so, following GL, we use 
the convective term to estimate the bulk kinetic dissipa- 
tion 



The dissipations estimated in Eons. 1251 and 1301 are 
sufficient to define the two relevant regimes for turbulent 
electroconvection. In one scenario, we assume that both 
the electrical and kinetic dissipations occur primarily in 
the boundary. Then Eon. I20| will be written as 



- 7 ^Re\ (25) 

The relations [21 and [53 are virtually identical to their 
counterparts for turbulent RBC given in Section 2.3 of 
Ref. 0. 



Electric dissipations 

The applied potential V drops over the electric poten- 
tial boundary length A^, which we assume is small such 
that <C ^i, <C To and <^ Tq ~ Vi. As with the 
kinetic dissipations, these requirements greatly simplify 
the averaging prescribed by Eans.l^andl^ Since most 
of the applied potential drops over A^ at the inner and 
outer electrodes, there is effectively an electrical short i.e. 
a constant potential, in the bulk. Invoking the standard 
arguments of RBC, where the thermal short determines 
a relation between the thermal boundary layer and the 
heat transport, we find that 



A,; 



{rj -|-ro)ln(l/a) 
Nu 



(26) 



Here, we have assumed that half the applied potential 
drops over the potential boundary layer at the inner elec- 
trode and half at the outer electrode. We assume that 
the electric potential boundary layer thicknesses are the 
same at the inner and outer electrodes. This symmetry 
is exact for radius ratio a ^ 1 or aspect ratio F — > cx). 
Our derivation implicitly assumes that the role of the 
asymmetry is not crucial at smaller F. 

In the boundary layer, the electric dissipation is given 

by 



(^(vv.)^)m^^ — 



A„ 



By balancing the latter two terms in Eqn. [3 



u • V?/' 



(27) 



(28) 



(29) 



we find, up to a linear expansion of gijjj), that 
V_ ^ cTiro-n)V 

By using U as the velocity scale, we have implicitly as- 
sumed that A^ < Au- Substituting in Eqn. [23 for V/X^ 
from Ean. 1291 and for A,^ from Eqn. [201 we get 



l + a\ ln(l/a) Re VaV^ 
1 — a J Nu {ro — ri)"^ 



(30) 



hi 



(31) 



This corresponds to the I; regime of Ref. 8]. For the 
left-hand-sides in the above equations, we use the glob- 
ally averaged dissipations derived from the equations of 
motion given in Eg ns . IT^ and l 1 91 For the right-hand-sides 
we use the boundary layer estimates given in Eons 1 2 41 and 
[2ni After some algebraic manipulation, we find 

Nu cx F{r)TZ^/^'P^/^ (regime I; of Ref. [8]). (32) 



Repeating the above procedure for 



bulk 



(33) 



which corresponds to regime 11; of Ref. j^l , we find that 

Nu oc F{T)'R}/^V^/^ (regime II; of Ref. [8]). (34) 

In general, this procedure gives a power law depen- 
dence for f in Eqn. [23 of the form 



f - TVV^ . 



(35) 



It also follows that the Reynolds number Re of the vor- 
tices, based on the radial velocity of the fluid and the 
separation between the electrodes — r^, scales as 

Rer^TV'V^'. (36) 

We find the same set of exponents 7,(5, 7*, (5* as the GL 
theory for turbulent RBC. The 1; and 11; regimes are 
relevant to the low TZ and broad V range that can be 
accessed by smectic electroconvection. 



aspect ratio dependence F(r) 

After balancing the dominant contributions, as in 
Eans.rmand l33l we find that the power law TVV^ has a 
common algebraic prefactor F that is only a function of 
the aspect ratio F. Unlike in previous studies of RBC in 
Cartesian geometry 38j, F(F) is not itself a power law. 
Instead, it is given by 



F(F) 




(37) 



F(F) specifies the aspect ratio dependence of the global 
charge or heat transport, as contained in the Nusselt 
number Nu. The Reynolds number Re of the large scale 
circulation, however, is local to each vortex and is inde- 
pendent of F. 
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EXPERIMENT AND RESULTS 

In this section, we describe the main features of the 
experimental apparatus and procedure, the data analysis 
and the results. Our apparatus is similar to one used 
previously for studies of thin film electroconvection in the 
weakly nonlinear regime 

I26I I271 1281 , which has more recently been adapted to the 
turbulent regime Other forms of electroconvection 
in bulk liquid crystals have similarly been extended to 
study scaling in the turbulent regime [13, 113| ■ 

The experiment consists of an annular liquid crystal 
film of octylcyanobiphenyl (8CB) freely suspended be- 
tween two concentric stainless steel electrodes as shown 
in Fig. ^ The annular film was about 2 cm in diameter. 
The experiment is enclosed by an aluminum box which 
serves both as a faraday cage and as a rough vacuum 
chamber. For the experiments discussed here, the film 
was at atmospheric pressure and was temperature con- 
trolled to 24 ± 2°C. At this temperature, 8CB is in the 
smectic A phase. In this phase, the elongated molecules 
align normal to the plane of the film, which consists of an 
integer number of layers. In 8CB, each layer is 3.16 nm 
thick [s^. In all experiments, we used films that were 
uniformly thick to within ±3 layers, and which had thick- 
nesses between 30 and 85 layers. Within the plane of the 
layers, the film closely approximates an ideal 2D Newto- 
nian fluid. The layered structure strongly restricts fluid 
motion perpendicular to the layers. The film thickness 
is comparable to the wavelength of visible light. We de- 
termined the thickness of a film by its interference color 
under reflected white light, using standard colorimetric 
functions [s^ [s^ . During the experiments we observed 
and recorded the reflected film color with a CCD camera. 

Pure 8CB has a low, uncontrolled electrical conduc- 
tivity due to residual ionic impurities. To control the 
conductivity, we dope the 8CB with tetracyanoquin- 
odimcthane (TCNQ) [s^l, an electron acceptor. With 
a concentration of TCNQ between 5 x 10~^ and 5 x 10^* 
by mass, we find that the sample has a bulk conductivity 
in the range 10^^ to 10^^ ft^^m^^. We determine the 
conductivity of each film from its ohmic response below 
the onset of convection, as discussed below. 

The experimental procedure consists of applying a DC 
voltage V across the film and measuring the resulting 
electrical current /. The inner (outer) electrode is elec- 
trically high (grounded). The applied voltage is incre- 
mented in small steps from to 1000 V and then decre- 
mented to zero. At each voltage, the current is mea- 
sured with a computer-interfaced electrometer, which is 
equipped with low noise triaxial cables. The film resis- 
tance is in the Tfi range and typical currents are ~ 1 pA. 
To determine the dimensionless charge transport Nu, we 
require an accurate value for the critical voltage Vc at 
the onset of convection. To measure Vc, we use a small 



voltage step, ~ 1 volt in the voltage range between and 
50 V, which brackets the typical critical voltage for most 
of our films. A larger voltage step of ^ 10 V is used in the 
range 60 to 1000 V. The larger step is necessary to limit 
the drift in electrical conductivity due to electrochemical 
reactions in the film. At each applied voltage, we make 
100 current measurements spaced by 25 ms. The average 
values of the current were used to calculate Nu, as de- 
scribed below. In addition to the IV data, we measured 
the film thickness s and aspect ratio F. Further details 
of the material preparation and experimental procedure 
can be found in Refs. [Ulillil. 

From the slope of the IV curve in the conduction 
regime where the film is quiescent, we determined the 
film conductance C. The Nusselt number is the current 
divided by the conduction current, A^w ~ I/Ic= I/CV. 
The control parameter TZ defined in Eqn. [SJ was calcu- 
lated using the accepted values of the viscosity, the mea- 
sured conductivity, film thickness and the applied volt- 
age. It varies between and ~ 10^. This is moderate 
compared with the very high Rayleigh numbers achiev- 
able in RBC. However, the critical value TZc of TZ at 
the onset of electroconvection is about a factor of 10 
smaller than that for RBC. The Prandtl number V de- 
fined in Eqn. was calculated from the material param- 
eters and film dimensions. In the experiments reported 
here, 5 < P < 250. 

Electroconvecting annular smectic films have several 
advantages for the study of turbulent scaling, relative to 
conventional RBC. In RBC experiments, the heat losses 
through the sidewalls must be taken into account and 
corrected for 0, whereas the annular film has no side- 
walls. The characteristic time scales of electroconvec- 
tion are many orders of magnitude shorter than in RBC, 
making data acquisition much faster. Also, the annu- 
lar aspect ratio can be easily varied over a broad range 
(0.3 < F < 17 for the experiments reported here). A 
similar range of F in RBC is possible in principle, but 
would be very cumbersome and time consuming in prac- 
tice. 

Electroconvecting smectic films also have some clear 
disadvantages, relative to RBC. Degradation of the liq- 
uid crystal under DC excitation results in significant con- 
ductivity drifts over an experimental run. The drift can 
be as large as 30% over the course of an experiment. The 
drift can be partially compensated for by monitoring the 
change in the critical voltage Vc before and after each 
sweep of the voltage, as discussed further below. The 
total DC charge flow is quite large; ~ 0.4/iC over a time 
interval of ~ Ih. The conductivity drift results in an 
uncontrolled variation of both TZ and T', increasing the 
uncertainty in these parameters. A second disadvantage 
is the rather modest upper limit on TZ, relative to that 
attainable in RBC. Higher TZ could be had by simply 
increasing the applied voltage V. However, large elec- 
tric fields will eventually lead to dielectric breakdown. 
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FIG. 3: A representative current-voltage curve for an annular 
film showing the onset of electroconvection. Data obtained 
for increasing (decreasing) voltages are shown in A(v)- The 
inset shows the rms fluctuations of the current vs. the applied 
vohage. Here T = 3.74 ± 0.02 and P = 23 ± 2. 

destroying the film. We estimate that with the present 
material parameters, the maximum accessible TZ is about 
10^. 

Fig.Oshows a representative IV response. The critical 
voltage Vc is identifiable by the upward kink in the IV 
curve at F ~ 40V, after which the film makes a transition 
from conduction to convection. When V < Vc, the fluid 
is motionless and charge is carried by ohmic conduction. 
When V > Vc, the electrical driving force overcomes dis- 
sipation and the film flows in a series of counter-rotating 
vortices, which are laminar for small V > Vc- The fluid 
circulation carries an additional current by convection, as 
is apparent by the increase in the slope of the IV curve 
above Vc- At higher voltages, the vortices become un- 
steady. The transition to unsteady flow is identified by 
a sudden jump in the current fluctuations at 200V, 
as shown in the inset of Fig. |3| With even higher driving 
voltage V > 600V, the flow becomes turbulent. It is in 
this highest range that Nu — TZ power law scalings were 
observed in the experiments. 

Dependence on the Rayleigh number TZ 

The electric Rayleigh number TZ, defined by Eqn.|Hl is 
a dimensionless measure of the external electric driving 
force. Stated in terms of experimentally measurable pa- 
rameters, TZ = ^oV^ / <73ri3s'^ , where eo is the permittivity 
of free space, V is the imposed voltage, 0-3 is the bulk 
conductivity, 773 is the bulk viscosity, and s is the film 
thickness. We find the bulk conductivity 0-3 by measuring 
the ohmic conductance of the film C = 2775(73/ In (l/a), 
where a is the radius ratio Vi/ro- C can be directly deter- 
mined from the IV curve, C = I/V, in the conductive 
regime for V < Vc- Previous experiments have deter- 
mined the bulk viscosity to be 773 = 0.19 ± 0.05 kg/ms at 



atmospheric pressure |25l |. 

As previously mentioned, we determine the film thick- 
ness s by matching the observed reflection color of the 
film to a color chart. 

The dimensionless Nusselt number Nu measures the 
convective contribution to the charge transport. Nu can 
be directly calculated from IV data; Nu is the total cur- 
rent normalized by the conductive current. The error in 
the scaling exponent 7 in the relation Nu ^ TZ'' stems 
from uncertainties in the film thickness s, the critical 
voltage Vc, and the film conductance C. The main source 
of error is the drift in the film conductance C. The drift 
results in slightly different slopes for the IV curves in the 
conductive regime between the increasing voltage and de- 
creasing voltage sweeps in a single experiment. The cause 
of the conductance drift is not well understood, but is 
presumably due to electrochemical changes in the liquid 
crystal material. A large drift in the film conductance can 
lead to a significant uncertainty in the Nu — TZ scaling. 
In principle, this drift could be completely compensated 
for if C were known for every IV measurement. Unfortu- 
nately, the conductance C cannot be independently de- 
termined from the data while the film is convecting. We 
must use the IV data in the conduction regime, which 
occurs before and after the convection regime during one 
voltage sweep. 

In order to bracket the drift, we have analyzed the 
data using two methods. From a fit to the ohmic re- 
sponse at the beginning and end of each sweep, we deter- 
mine two conductances. Cup and Cdown- In method A, 
we use Cup to reduce the data obtained when increment- 
ing the voltage and Cdown to reduce the data obtained 
when decrementing the voltage. This method concen- 
trates the conductance error at the maximum voltage, 
and thus overestimates the effect of drift at the highest 
voltages. In the second method, which we refer to below 
as method B, we instead assume that the conductance C 
varies linearly between Cup and Cdown while the fluid is 
convecting. This method can be thought of as a linear 
approximation to the unknown evolution of C while the 
film is convecting. It probably underestimates the true 
drift. By comparing our results using methods A and 
B, we can gauge the overall effect of the drift. We have 
analyzed all the data using both methods. In some runs, 
especially large drifts made compensating for the drift 
essentially impossible. We discarded runs in which the 
difference in Cup and Cdown exceeded 30%. 

Fig. 0Jl shows Nu vs. TZ data for three aspect ratios 
r. When Nu — 1, the fluid is quiescent. When Nu > 1, 
the fluid is convecting. For TZ > 10*, the data reveal 
power law behavior with Nu ^ TZ'^ and values of 7 close 
to either 1/5 (for smaller V) or 1/4 (for larger V), in a 
good agreement with the theoretical predictions for the 
small TZ regime. Detailed experimental results for various 
F and 'P are listed in Tabled The Rayleigh number TZ 
range used to fit Nu ~ TZ'^ was between approximately 
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TABLE L Results of fits to iV-u ~ 7^^, for different aspect 
ratios F and various Prandtl numbers V by analysis method 
A. Results obtained by method B are consistent. 



r 


range of T' 




range 


of 7^ 




0.33 


6-9 


0.21 ± 0.02 


8 X 10^ 


- 7 X 


10* 


1.54 


19-28 


0.21 ± 0.01 


1 X 10* 


- 2 X 


10^ 


3.74 


21 - 25 


0.22 ±0.02 


2 X 10* 


- 2 X 


10^ 


6.6 


25 - 41 


0.24 ± 0.04 


2 X 10* 


- 2 X 


10'^ 


6.6 


48-61 


0.21 ±0.01 


2 X 10* 


- 3 X 


10'' 


11.1 


70-74 


0.25 ± 0.02 


1 X 10* 


- 2 X 


10^ 


11.1 


112 - 120 


0.18 ±0.01 


2 X 10* 


- 4 X 


10'^ 


11.1 


127 - 136 


0.19 ±0.01 


2 X 10* 


- 4 X 


10^ 


16.1 


205 - 241 


0.18 ±0.02 


8 X 10^ 


- 5 X 


10'^ 



100 X TZc and the final data point at the largest voh- 
age, 1000 volts. Also note that for annular electrocon- 
vection, the critical Rayleigh number TZc at the onset of 
convection is ^ 100, and thus a factor of 10 smaller than 
the corresponding Rac ^ 2000 for RBC. From a total 
of 46 experiments, the scaling exponents 7 were either 
0.20 ± 0.03 (0.19 ± 0.03) or 0.25 ± 0.02 (0.24 ± 0.02) by 
analysis method A (B), depending on V. The error bars 
quoted cover the scatter in 7 and also include the uncer- 
tainties in the film thickness, the critical voltage, and the 
film conductance drift. In Fig. compensated plots of 
Nu/7U vs. TZ are shown. The best fit exponent 7 ob- 
tained from experiments was used to compensate the Nu 
data. Although only about one decade of scaling range is 
available, the compensated plots show that a local power 
law is an adequate description of the data. This range 
is probably too short to resolve details about possible 
cross-over scalings ^] . 



Dependence on the Prandtl number V 

The electric Prandtl number V defined by Eqn. is 
the dimensionless ratio of the charge and viscous relax- 
ation time scales. V^^ appears as a prefactor in the 
nonlinear and time derivative terms when the equations 
of motion are written in dimensionless variables. It is 
thus reasonable that any dependence on V vanishes for 
large V. However, for P ~ 1 the turbulent flow and 
scalings are expected to depend on V. The relative 
length scales of the electric potential and viscous bound- 
ary layers, which depend on enter into the scaling 
arguments of Grossmann and Lohse 0|. For RBC, GL 
theory predicts that Nu should exhibit local power law 
scalings with Prandtl number, albeit with rather small 
powers 0, 13 . RBC experiments suggest that the heat 
transport is largely independent of the Prandtl number 
between 4 to 1350. For example, it has been found that 
Nu ^ -p-o-os at r = 1 |1J]. In turbulent electrocon- 




FIG. 4: (Color online) (a) shows representative plots of Nu vs 
7^ for r = 0.33 ±0.01 (A), 3.74 ±0.02 (□), and 6.60±0.02 (o), 
analyzed by method A. For these, V = 8.8 ± 0.5, 21 ± 1, 
and 36 ± 1, respectively. Least-square fits to the power law 
Nu - give best fit values 7 = 0.19 ± 0.01, 0.21 ± 0.02, 
and 0.26 ± 0.01, respectively. The solid reference lines in the 
figure have slopes of 1/5 (lower one) and 1/4 (upper one). In 
(b) we plot the compensated scaling Nu TZ ' vs. TZ for the 
same data as in (a) . The inset in (b) shows a more expanded 
scale for F = 0.33 (A), and 3.74 (□). 

vection for 2 x 10'* <TZ< 10^, we find that Nu varies by 
only a factor of 2 over the broad range 5 <P < 250. 

Fig. O shows a plot of Nu vs. V for two aspect ratios 
r = 6.6 and 11.1. At each V, we have averaged the 
compensated data NuTZ'^ over 2 x 10" < 7^ < 10^, with 
7 from the best fit. The data plotted in Fig. |31 were 
obtained from a total of 26 experiments, and results for 
both analysis methods A and B are shown. 

For a fixed aspect ratio F, our data suggest a crossover 
from one local power law Nu ~ TZ"' to another as V 
increases. Taking T = 6.6 as an example, the Nu vs. 
TZ scaling for 20 < "P < 40 gives a 7 exponent ~ 1/4 , 
while for 50 < < 70, the exponent is « 1/5. The same 
indications of a 1/4 to 1/5 crossover of scaling exponents 
are found for F = 11.1, but for a higher value of "P. In 
the case of F = 11.1, we find 7 « 1/4 for 70 < P < 80 
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and 7 « 1/5 for 100 < V < 130. One interpretation of 
this observation is that the boundaries between different 
Nu vs. TZ scaUng regimes depend on the aspect ratio 
r. This has not been previously considered in scahng 
theories of RBC, which are specific to the case F ~ 1. 
In annular electroconvection, we may have to consider a 
three dimensional parameter space of Nu scaling regimes 
that depends on TZ, V and F. 

The GL scaling theory applied to annular electrocon- 
vection yielded an explicit aspect ratio dependence for 
the dimensionless charge transport Nu, given by the 
function F(F) defined in Eqn. |^ Thus, we can take 
the varying aspect ratio into account by dividing the Nu 
data by the theoretically predicted value of i^(r). FiglSJ 
shows such a fully compensated plot of NuTZ 'F^^ vs. 
V for various aspect ratios F. Our data span 5 < P < 250 
and 0.3 < F < 17. All of this data had 7 ~ 1/5. A weak 
dependence on V remains, amid considerable scatter. A 
power law fit - V^^ gives (3 = 0.20ib0.04 (analysis method 
A) and 0.26 ± 0.05 (analysis method B). The theoretical 
prediction is Nu oc F{T)TZ^^^V^^^ for the regime of low 
TZ and Nu. Thus, our measured value of P is at least 
consistent with the GL theory applied to annular elec- 
troconvection for this specific regime, although the scat- 
ter is obviously too large to definitively establish that a 
power law is present. For the the adjacent regime where 
J — 1/4, we have less than a decade of T^ range which is 
insufficient even to look for such consistency. 

Dependence on the aspect ratio F 

In the scaling theory described above, we explicitly 
accounted for the aspect ratio dependence, a considera- 
tion omitted in the GL theory for turbulent RBC 0, ^ 
which treats the case F = 1 . In our formulation of the 
GL theory for annular electroconvection, we found that 
the charge transport is modified by an aspect ratio de- 
pendent prefactor F(r) given by Eqn. |23 Unlike the 
previous studies of RBC in Ref. [s^, we find the aspect 
ratio dependence is not a power law scaling but rather a 
simple function of the annular geometry. To make a di- 
rect comparison to previous turbulent RBC experiments, 
we consider a new function fcF, which we define to be 
F(F) multiplied by a normalization constant k, chosen 
such that fcF(F = 1) = 1. The appropriate value of k 
is 7r/((7r + l)ln(27r + 1)) = 0.382. The function A:F(F) 
decreases monotonically with F with its greatest varia- 
tion for F < 2, and is within 2% of its limiting value 
fcF(oo) = 0.764 for F > 7. 

Our experimental data span the range 0.3 < F < 17. 
Because our data span the wide range 5 < P < 250, we 
expect some corrections due to changes in the Nu vs. T' 
scaling, predicted by the theory and also observed in the 
current experiments. To separate the aspect ratio de- 
pendence of Nu and to compare it with the theoretical 
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FIG. 5: (Color online) Plots of averaged Nu/Tf vs. V for 
r = 6.60±0.05 (black symbols) and ll.liO.l (blue symbols), 
where 7 is taken from the best fit to Nu ~ TV . Circular 
symbols (o, •) show results obtained by analysis method A, 
while (*) symbols are obtained by method B. Solid (•) and 
open (o) symbols indicate when the scaling exponent 7 ~ 
1/4 (1/5). Note that the crossover between 1/4 and 1/5 
exponents occurs at different V for different aspect ratios F, 
as indicated by the shaded rectangles. 
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FIG. 6; (Color online) The average value of fully compen- 
sated Ai'it/(F(r)7?.^) vs. P on a log-log scale for various 
r. Circular symbols (o) show results obtained by analysis 
method A, while (*) symbols are obtained by method B. All of 
these data had a Nu vs. TZ scaling exponent of 7 ~ 1/5. The 
solid (dashed) line is the best fit to ~ P'' with (5 = 0.20 (0.26), 
using analysis method A (B). 



prediction Nu oc F(F)7^''7"^, we divide Nu by TZ''T'^ and 
take its averaged value over the TZ range from 2 x 10^ 
to 10^. Nu could still depend on F independent prefac- 
tors which are not captured by a scaling theory. These 
prefactors can not be separately extracted from the ex- 
perimental data and may change for the various scaling 
regimes jsj. To avoid these, we restrict our discussion 
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to those data for which the Nu vs. TZ scahng exponent 
7 is close to 1/5. We extract the exponents 7 from the 
power law fits of Nu ^ TZ"^ to the experimental data. We 
use 5 = 0.20 and 0.26 which are obtained by the power 
law fit of 7V■u/(F(^)7^T) - shown in Fig. using 
the two analysis methods A and B, as described above. 
This effectively completes a circle of mutual consistency 
checks which delivers a self-consistent experimental result 
for the F(r) dependence alone. Using one free parameter 
for all the data, we scale these data so that F(r — 1) — 1, 
to allow a comparison to similarly normalized RBC data. 
Our data for six different F, obtained from a total of 36 
experiments, are in a reasonable agreement with the the- 
oretical prediction for fcF, as shown in Fig. [3 The errors 
are representative of the scatter in Nu /TW^ among runs 
for each aspect ratio. 

The data at F = 6.6 and 16.1, which deviate most 
from the theoretical prediction for kF, consist of a few 
runs which have rather high values of V, about 60 and 
220, respectively. This deviation may be due to difhculty 
of finding the appropriate V scaling exponent 5. For rel- 
atively small TZ and high V, scaling theory predicts 
that Nu becomes independent of V, scaling with TZ^/^V^. 
It is possible that our few data points with very high V 
might fall into the scaling regime 7 = 1/5, S — instead 
of 7 = 1/5, 5 = 1/5, as we have assumed. Thus, the 
F(r) data at F = 6.6 and 16.1 may be significantly un- 
derestimated in Fig. \7\ Many more experiments with a 
wider range of parameters would be required to system- 
atically explore the various scaling regimes, boundaries 
and crossover effects. 

Data from several turbulent RBC experiments [E^IE 
[15I Is^ for various values of F are also in broad agree- 
ment with the function fcF(F), in spite of the difference 
in geometry and the higher range of Rayleigh numbers. 
Details of this comparison may be found in Ref. [f^ . 
One should only expect the comparison between the as- 
pect ratio dependence for RBC and the function kF to 
be reasonable in the limit F — > 00. Nevertheless, we find 
approximate agreement, in spite of the difference in ge- 
ometry, the smaller aspect ratios and the higher range of 
Rayleigh numbers. 

CONCLUSION 

We have theoretically and experimentally studied how 
the dimensionless charge transport Nu scales with the 
Rayleigh number TZ, the Prandtl number and the as- 
pect ratio F, for turbulent electroconvection in a 2D an- 
nular film. The electroconvection is driven by an unsta- 
ble charge distribution that is analogous to the inverted 
fluid density distribution in 3D, thermally driven RBC. 
The strong similarity of the governing equations between 
electroconvection and RBC allowed us to adapt GL scal- 
ing theory |^ and experimentally investigate its con- 
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FIG. 7: Plot of fcF(r) vs. T from the scaling theory (solid 
line), with k chosen so that fcF(r = 1) is unity. Experimental 
data for turbulent electroconvection are shown by solid sym- 
bols (•) by method A and by (*) for method B. The data span 
5 < P < 250. The inset shows the same data on a logarithmic 
scale. 

sequences. The unique annular geometry of the electro- 
convection also made it possible to explicitly account for 
the aspect ratio dependence of the scaling relations. 

From the theory, we found various regime-dependent 
local power laws of the form Nu ~ F{T)TZ''T^^ , with the 
same exponents as those for turbulent RBC, but with an 
additional aspect ratio dependent prefactor F(F). F(F) 
is a nontrivial function of the finite annular geometry, 
rather than a power law. 

In experiments, we found that the exponents for the 
Nu vs. TZ scaling were consistent with 1/4 or 1/5, for 
10** ^TZ <2x 10^, in reasonable agreement with the the- 
ory for the regime where TZ and Nu are both small. The 
Nu measurements for 0.3 < F < 17 and 5 < < 250 are 
consistent with the theoretical prediction for the regime 
where Nu F{T)TZ^^^T^^^^ . Furthermore, our exper- 
imental data suggests F-dependent cross-over between 
different scaling regimes. 

The weak dependence of F on F for large F suggests 
that the global heat or charge transport approaches a 
universal, F independent limit for laterally extended sys- 
tems. This conjecture suggests that future work on the 
charge or heat transport in turbulent convection should 
focus on large F systems. Although such systems can be 
difficult to achieve experimentally, they are common in 
nature. 

The parameter space of TZ, "P and F is, however, very 
large. In this study, we have have only sparingly sam- 
pled from the experimentally accessible portion of this 
parameter space. In spite of the general consistency be- 
tween the experiment and theory, it is difficult to draw 
definitive conclusions about the scaling assumptions un- 
derlying the theory. The range of parameters we can 
access experimentally is too narrow to span the various 
regimes. Also, the decomposition of the dissipations into 
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bulk and boundary contributions is not directly testable 
with our current experimental techniques. 

The great strength of the GL scaling analysis is its 
generality. It can be applied as easily to our 2D elec- 
troconvection as 3D RBC. The relative simplicity of the 
2D fluid mechanics in electroconvection suggest that we 
can greatly extend the range of parameters by numerical 
simulation. Simulations will also allow us to test directly 
the scaling assumptions. The GL theory could also be 
extended to the case of convection with a superposed 
shear [2^ 123 . I27I , a situation which is experimen- 
tally feasible in electroconvection but not in RBC. 

The nearly power law scaling of the globally averaged 
heat or charge transport in turbulent convection poses 
an interesting and difficult problem. Any approach to 
this problem must creatively combine theory, simulation 
and precision experiment. The unique features of thin 
film electroconvection give us a new vantage point on 
this challenging, and as yet unsolved problem. 

This research was supported by the Natural Science 
and Engineering Research Council of Canada. 
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